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Abstract 

We discuss system with spin 5 = 3/2 with general isotropic nearest- 
neighbor exchange within a mean-field approximation possess. we drive 
equations describing non-Heisenberg isotropic model using coherent states 
of SU(4) group in real parameters and then obtain dispersion equations of 
spin wave of dipole and quadrupole branches for a small linear excitation 
from the ground state. 

1 Introduction 

Magnetically ordered materials (magnets) are known as essentially nonlinear 
systems [1,2] . Localized nonlinear excitations with finite energy, or solitons, 
play an important role in description of nonlinear dynamics, in particular, spin 
dynamics for low-dimensional magnets, with different kind of magnetic order. 
To date, solitons in Heisenberg ferromagnets, whose macroscopic dynamics are 
described by the Landau-Lifshitz equation for the constant-length magnetization 
vector, have been studied in details. [2,3,4,5,6]. In terms of microscopic spin 
models, this picture corresponds to the exchange Heisenberg Hamiltonian, with 
the isotropic bilinear spin interaction J (5*1.52), where J is exchange constant 
and i,j are nearest-neighbor in corresponding site . For a spin of 5 > 1/2, 
specially, for Single-Ion exchange, for more carefully description of magnets, the 
isotropic interaction is not limited by this term and can include higher invariants 
such as J(Sx.S 2 ) n with n up to 2S [7]. 

Given that each spin state can be written in the form — X^i" 1 IV**)) then 
25 + 1 complex parameters are required to describe each spin state and it is 
equivalent to 45 + 2 degrees of freedom . However, one degree for normalization 
condition, (ip\ip) = 1 , and another one for arbitrary phase decrease. At the 
end, 45 degrees of freedom are required to describe spin state fully. [8] 
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In this paper, the general isotropic model with the spin S = 3/2 and the 
nearest neighbor interaction is described by the Hamiltonian 

H = - J2(J(Si.S i+ i) + K(Si.S i+1 ) 2 + L(£4 +1 ) 3 (1) 

i 

Where Si are spin-3/2 operators at the lattice site i ; J, K, and L are the 
exchange constants, corresponding to the bilinear, biquadratic, and the bicubic 
exchange interactions, respectively, and summation over pares of the nearest 
neighbors is implied. This Hamiltonian defines a one-dimensional ferromag- 
netic spin chain and J is a positive coefficient, also, coefficient K is exchange 
integral for quadrapole moment and coefficient L is exchange integral for oc- 
tapole moment. 

The main purpose here is to develop classical equations for Hamiltonian (1) 
and then finding solutions of spin wave for small linear excitations above the 
ground state. We know that coherent states minimize the uncertainty relation, 
then these states arc the nearest states to classical pictures. To this end, in 
section 2, coherent states for spin S — 3/2 which are the coherent states in 
SU(4) group are introduced. To obtain classical Hamiltonian, average values of 
spin operator are needed, so in section 3, these values and classical spin relation 
are derived. In section 4, the resulted Hamiltonian is substituted in classical 
equation of motion obtained by acting Feynman path integral over coherent 
states. Finally, obtaining spin wave equations and dispersion equations in case of 
small linear excitation from the ground state for dipole and quadrapole branches 
put an end to this research. 

2 Coherent states for spin 5 = 3/2 

Coherent states are special kind of quantum states that their dynamics are very 
similar to their corresponding classical system. That which coherent state can be 
used in problems depends on the operators symmetry in those peoblcms. Due to 
the operators symmetry in Hamiltonian (1), for full description and considering 
all multipole excitation, we used from coherent states in SU(4) group, in this 
group vacuum state is (I,0,0,0) T and its coherent state is introduced as: [9] 

= D 3 / 2 (^,7)e 2 ^V^V ife# * V >> 

(2) 

In above relation D 3 / 2 (#, <p, 7) is Wigner function for spin S = 3/2, operator 
Q xy is related to quadrupole moment and operator F xyz is relevant to octapole 
moment. 

This six-parameters state have the properties of SU(4) coherent state. Two 
angle, 9 and <j>, Euler angles, define the orientation of the classical spin vector. 
The angle 7 is the rotation of the quadruple moment about the spin vector. 
The parameter, g, defines change of the spin vector magnitude and that of the 
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quadruple moment, the angle /3 is the rotation of the octuple moment about the 
spin vector and parameter k, defines change of the spin vector magnitude and 
that of the octuple moment. It should be noted that, range of angles, <f), 7 and 
(3 are between zero to 2ir and angle 9 change between -7r to it. Using the state 
(2), one can construct the coherent state path integral and find the lagrangian 
L:[9] 

L = cos2kcos 2 g(3cos 2 g(3 t + cos6<j) t + 7*) - H(9, <j>, g, 7, /3, k) (3) 

Where x t — d/dt and H is the classical (mean- field) energy of the system, 
which equals to the quantum mean value of the Hamiltonian with the state (2). 
Note that in deriving Lagrangian of spin system using path integral, two more 
terms appeared. First term is dynamic term that has Berry phase characteristics 
and takes many attentions in phenomena like spin tunneling and the other term 
is boundary term depends on boundary values of path. Both of these terms are 
not interested here and hence omitted. 

3 Average spin operators and their products in 
SU(4) group 

Here we consider classical counterparts of the spin operators and their products 
contained in the Hamiltonian (1) . The vector 

s=&0m (4) 

Can be regarded as a classical spin vector, and 

Q« = I(54+^-^7) (5) 

components of quadrupole moment. Because we can write any coherent state 
as multiple of single site coherent states, namely: 

Then Spin operators in ground state of non-single ions Hamiltonian can be 
commute in different lattices [10]; so 

M&AM) = M& n \m\Si +1 W) (7) 
where \ip) = n \ip) n+1 . 

Then average values expressions in SU(4) group for coherent state (2) are: 
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(5+) = |e^(l- Acos 2 g)cos2ksin6 

(S~) = |e"^(l- 4:COS 2 g)cos2ksin6 

(S z ) = ^(l-4cos 2 g)cos2kcos0 (8) 



and also 



S 2 = ^(l-4 COS 2 . g )W2fc 
9 

g 2 = -(1 - 4cos 2 g) 2 sin 2 2k 

S 2 + q 2 = \{l-Aco S 2 g) 2 (9) 

In the above relation S 2 related to dipole moment and q 2 is related to 
quadrapole moment, if we only cosider quadruple moment, must be set 5 = 0, 
then 



S 2 + q A = cte (10) 

In classical limit, if we used from above relations for Hamiltonian (1), the 
classical Hamiltonian in continous limit obtained in the following form: 

= _ f ^^jA 2 co S 2 2k+^KA i cos i 2k+ 7 -^-LA 6 cos 6 2k 
J a M 16 64 

+ -KA 2 cos 2 2k + —LA 4 cos 4 2k) x 
2 y 2 16 ; 

9 

x(36glcos 2 2ksin 2 2g + -A 2 cos 2 2k(6 2 x + sin 2 Bc\> 2 x ) 
+§k 2 x A 2 sin 2 2k)) (11) 

where ao is a lattice crystal length and 



A = (1 - 4cos 2 g) (12) 

4 Classical Equations 

To obtain the dynamic equation, we vary the Lagrangian (3), obtaining the set of 
equations, but because we only consider quadruple moment we set 5 = 0. These 
equations completely describe nonlinear dynamics of Hamiltonian of problem 
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under study up to quadrupole branch. Solution ol these equations are magnetic 
solitons. then classical equation is: 



1 81 

— 4> t = (cos2fc(128J + 81(32if + 729L) + 32A(8K + 243L)cos4fc 

luq 512 

xl9683Lcos8k)(6 xx csc6 + (f> 2 x cos9)al 
1 81 

— 6 t = (f> xx cos2k(16J + 6A8Kcos 2 2k + 19683Lcos 4 2fc) 

wo 64 

xsinOciQ 

1 27 

— p t = — (16 Jcos2k + 6A8Kcos 3 2k + 19683Lcos 5 2fc) 
wo 32 

27 

H (39366A; xx cos 4 2fcsm2/c + 8k xx sin2k(AJ + 81K + 8\KcosAk) 

64 

+19683icos 5 2fc(4fc2 - 30* - 3<j) 2 x sin 2 6) + 6A8cos 3 2k 
x{AKk 2 x - 2K6 2 X - 2A3Lk 2 x sin 2 2k - 2K (gsin 2 6) - 16cos2fc 
x {\Q2Kk 2 x sin 2 2k + J{-Ak 2 x + 9 2 x + ^ 2 x sin 2 9)))al 
1 81 

— 7t = (16 Jcos2A; + 648Kcos 3 2fc + 19683icos 5 2fc) 

w 16 

81 1 

- — (19683Lcos 5 2A;(8A;2 - W 2 X + -(-5 + 7cos20 + 6 xx cot6)4> x ) 

+6AJk xx sin2k + 78732Lk xx sin2kcos 5 2k + 8cos2k(J (l(Sk 2 x - A0 2 x 
-<j>l + 3(f> 2 x cos29 + 29 xx ) - 6A8Kk 2 x sin 2 2k) + 6A8cos 3 2k 

{8Kk 2 x - AK6 2 X + K(^<f>l(-3 + 5co.s26>) + 6 xx cot0) - A8&Lk 2 x sin 2 2k) 
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+6A8Kk xx (sin2k + sin6k))al 
■kt = 



OJ 

—9t = (13) 

LJQ 

Where uj n — hao- In the above equations, if quadrapole and octapole excita- 
tions are ignored, (g = 0, k = 0), these equations are redused to Landau-Lifshitz 
equation. As a result, these equations compared to Landau-Lifshitz equation 
are more complete and also they enjoy higher degrees of freedom. 

In order to investigate the ferromagnets with anisotropy, it is necessary to 
find the classical ground states of these magnets. In this order, we consider in 
Hamiltonian (11) only part of without derivative: 

/At 9 9 9 

— {Q 2 Jcos 2 2k + Q 4 Kcos 4 2k + Q 6 Lcos 6 2k) (14) 
ao 2 2 2 

In order to this terms is minimum must be k = or k = ^ . The minimum 
value of classical Hamiltonian (14) is 
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ho = --((h 2 J+& i K+( 9 -r L ) (15) 

CLq z z z 

Let us examine the classical equations near the ground state, in this order 
we change the following relation for the paremeter k 



2fc^7T + fc (16) 

in the point 6 = | , in linear small excitation the classical equations change 
in the following form 

—<j) t = A9 xx al 

—Ot = -A(f) xx al 

UJQ 

1 8 2A 

Pt - -7T 

uj 3 

—It = -4A 

—h = o 

uj 

1 -9t = (17) 



OJQ 



where 



A = — (16J + 648X+ 19683L) (18) 

and x t = ^. Although exchange integrals K and L are small, numerical val- 
ues which are multiplied by them in the above relation demounstrate that 
quadrapole and octapole excitations in this problem are noticeable. Also, it 
is shown that for isotropic ferromagnets, the magnitude of average quadrupole 
moment is constant and its dynamics, is rotational dynamics around the classi- 
cal spin vector . 

We consider now the dispersion of spin wave propagating near the ground 
state. To this end we considering two functions 0, <f> in the form of plane waves, 

<t> = (t) e i{u)t ~ kx) + (t>Qe' l(ult - kx) 

= 6 Q e i(ut - kx) + 6 e- i{ut - kx) (19) 

we obtain the following equation for the spin wave propagating near the 
ground state for this isotropic Hamiltonian 
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= (J + 2K + 'iLfalulk* 
2 . 

uj 2 = -Aw 

lu 3 = -4Alo (20) 

It is evident from equation (20) that in addition to the dispersion acoustic 
branch, there exist two non-dispersion optical branches which are related to the 
dipole and quadrapole excitations. 

5 Discussion 

In terms of spin coherent states we have investigate S = 3/2 spin quantum 
system with the bilinear, biquadratic and bicubic isotropic exchange in the 
continuum limit, the proper Hamiltonian of the model can be written as bilinear 
on the generator SU(4) group[ll]. Knowledge of such group structure enables 
us to obtain some new exact analytical results. The analysis of the proper 
classical model allows us to get different soliton solutions with hnite energy and 
the spatial distribution of spin-dipole, spin-quadrupole and/or spin-octapole 
moments termed as dipole, quadrupole, octapole and dipole-quadrupole, dipole- 
octapole and quadrupole-octapole soliton, respectively. 
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